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SUMARY 

This  paper  investigates  mathematically  the  effect  on  the  pitch-yaw 
oscillations  of  a  spinning  re-entry  body,  which  will  arise  as  a  consequence 
of  the  body  having  a  small  amount  of  asymmetry. 

It  shows  how  this  effect  can  be  related  to  an  equivalent  re-entry 
disturbance  for  a  corresponding  symmetrical  body. 

Some  numerical  checks,  obtained  on  a  digital  calculating  machine,  are 
included,  as  substantiation  of  the  mathematical  theory. 
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1  INTRODUCTION 


The  oscillatory  motion  of  a  spinning  body,  descending  through  the  earth's 
atmosphere,  has  been  the  subject  of  a  fair  amount  of  theoretical  research  over 
the  last  few  years.  The  problem  is  to  derive  expressions,  preferably  reason¬ 
ably  simple  ones,  which  will  define  the  oscillatory  motion  in  terms  of  two 
groups  of  parameters. 

(1)  the  body's  inherent  characteristics  -  mass,  moments  of  inertia, 
aerodynamic  derivatives,  etc., 

(2)  the  flight  conditions  at  re-entry  -  velocity,  spin  rate,  the 
body's  orientation,  etc. 

This  is  an  important  issue,  because  the  oscillatory  motion  of  a  re-entry 
body  has  various  practical  repercussions.  It  governs,  for  example,  the  magni¬ 
tude  of  the  lateral  loading  experienced  by  the  body  during  descent  and  this  may 
be  a  critical  design  feature.  It  is  also  essential  to  have  a  theoretical 
insight  into  this  problem  in  order  to  evaluate  flight  data  (e.g.  accelerometer 
recordings)  obtained  during  a  re-entry  trial,  as  for  example  with  Black  Knight. 

In  Ref.1,  the  re-entry  oscillation  problem  is  treated  more  comprehensively, 
so  it  is  believed,  than  elsewhere.  There  are  inevitably,  however,  various 
additions  and  refinements  which  might  in  principle  be  made.  One  in  particular 
which  has  occasionally  given  rise  to  speculation,  is  what  would  happen  if  the 
re-entry  body  were  not  quite  symmetrical.  Ref.1  assumes  a  perfectly  symmetrical 
body  which  is,  of  course,  normally  the  designer's  intention,  but,  in  praotice, 
there  may  well  be  some  departure  from  this.  The  present  paper  deals  with  the 
problem  of  calculating  the  effect  on  the  oscillatory  motion  which  will  result 
from  such  asymmetry. 

It  i3  not  difficult  to  see  that  this  might  well  be  important.  Any 
asymmetry  in  a  spinning  body  will  constitute  a  forcing  input  at  the  spin  fre¬ 
quency,  this  frequency  being  substantially  constant  throughout.  The  natural 
frequency  of  the  body's  oscillations,  however,  increases  steadily  during  re¬ 
entry  (until  a  time  shortly  before  ground  level)  due  to  the  growth  of  the 
dynamic  pressure,  'g’pV^,  in  the  aerodynamic  forces.  Almost  always  this  natural 
frequency  varies  during  descent  from  well  below  the  spin  rate  to  well  above  it. 
The  system  passes,  therefore,  through  a  state  of  resonance.  As  is  well  known, 
if  a  forcing  condition  at  or  near  resonance  of  an  oscillating  system  is  main¬ 
tained,  large  amplitudes  of  motion  will  ensue.  The  general  issue  which  arises 
in  the  present  problem  is  to  discover  how  the  response  of  a  forced  oscillatory 
system  is  modified  by  the  continuous  state  of  flux  of  some  of  the  system's 
parameters. 

2  MATHEMATICAL  DEFINITION  OF  THE  PROBLEM 


Asymmetry  in  a  re-entry  vehicle  can  take  different  forms  (c.f.  Ref. 2). 
We  are  concerned  here  with  aerodynamic  asymmetry,  that  is  to  say  a  case  in 
which  the  roll  inertial  axis  is  inclined  at  a  small  angle  (e)  to  the 
aerodynamic  axis.  The  moments  of  inertia  in  pitch  and  yaw  are  equal. 
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We  follow  the  same  approach  as  in  Ref.1,  so  that,  as  seen  from  in  front 
with  small  perturbations  we  have  the  situation: - 


r 


a  aerodynamic  axis 

i  inertial  axis 

a,p  components  of  incidence  of  inertial  axis 

Y,p  angular  perturbation  components  of  the  trajectory 
q,r  angular  rate  components  of  inertial  axis 
e  angle  between  a  and  i 

A  roll  moment  of  inertia 

B  pitch  (or  yaw)  moment  of  inertia 

m  mass 

P  spin  rate 

The  equations  of  angular  motion  of  the  inertial  axis  are:- 


m  |_q+  ir  -  a-  ip 


=  ipVSC 


r~ 

I  a-  ip  + 

'al¬ 


ee 


iPt 


(1) 


— »- _  -  2 

Bq+ir-iAPq+ir  =  ■gpV  S£ 


iPt 


C  f  a-  ip  +  ee  +  C  (  q+ ir  +  i 

L  ma  V  )  mq  \ 


.  „  iPt' 
lePe 

...  (2) 
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These  equations  correspond  to  (7)  and  (8)  in  Ref.1,  with  the  additional  terms 
in  the  aerodynamic  forces  arising  from  the  asymmetry.  We  now  eliminate  q+  ir, 
to  obtain  a  differential  equation  in  a  -  ip. 


As  is  brought  out  in  Ref.1,  various  terms  appear  which  are  quite 
negligible  in  practice.  Retaining  only  terms  of  numerioal  significance,  we 
obtain: - 


a-  ip  +  a  -  ip 


.  AP 

iT 


+  a-  ip 


^pV2S€C 


m 


a 


B 


^PvscL 


m 


ipV2S£C 


m 


B 


a  iPt 
ee 


(3) 


This  is  the  counterpart  of  equation  (12)  in  Ref. 1,  the  asymmetry  appearing  on 
the  right-hand  side  of  the  equation. 

During  a  re-entry  descent,  the  velocity  is  substantially  constant  until 
the  lower  regions  of  the  atmosphere  are  reached,  whilst  the  air  density  builds 
up  exponentially  with  loss  of  height.  For  much  of  the  desoent,  therefore, 
p  takes  the  form 


=  p  e 
ro 


Kt 


(pQ,  K  are  constants) 


Restricting  the  investigation  to  this  case,  in  equation  (3)  we  replace 

v,  Kt 

p  by  pQ  e 


CO 


■4pV2S£C. 


ma  ,  2  Kt 

—  by  w  e 


CT  1  Cm 

L  m 

a  =  -kiVS  (  — - -~2- 

*-r  V  m 


B 


by  a  e 
o 


Kt 


CT  l  C 
/  m 

tnVR  ! -  x  - —3. 


b  =  ipvs  — -  +  B 


.  .  Kt 

by  b  e  • 
o 


We  then  arrive  at  the  differential  equation 
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a  -  ij3+  a-  ip 


Kt  .  AP 

V  '  1 T 


+  a-  ip 


2  Kt  .  AP  ,  .  ^Kt 

oo  e  -  i  —  (a  +  b  )e 
o  2B  '  o  o' 


2  Kt  iPt 
-eoo  e  e  . 
o 


This  is  the  equation  whose  solution  is  required. 
3  MATHEMATICAL  SOLUTION  OF  THE  PROBLEM 
We  make  a  change  of  variables 


...  (4) 


Kt 

T  =  e2  ,  a-  ip  =  Tiv  e 


a  T 
o 

2K 


,  AP 

T)  ,  where  v  =  — 


and  we  write 


Then  equation  (4)  transforms,  mathematically  into 


dT2  41  Lk2 


f4!  211>"  APP 

K2  '  K2  B 


-V  - 

K 


Tl 


=  -  C  ^  j2+i(2X-v)  e  2K 
K2 


(5) 


3.1  Special  case  of  no  damping  or  lift 


It  will  simplify  the  mathematical  development  of  the  solution  of  the 
general  equation  if  we  consider  first  the  special  case  in  which  there  is 
neither  aerodynamic  damping  nor  lift,  namely  the  case  in  which  the  aerodynamic 
forces  on  the  body  consist  solely  of  a  restoring  moment.  In  this  case,  aQ  and 

bQ  are  both  zero  and  equation  (5)  becomes 


T 


d2n 

dT2 


dT 


_  K 


T2+i(2X-v) 


(6) 


Now  make  the  further  substitution 
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ao  T 


o 


K 


z 


in  which  case  the  equation  becomes 


dz 


(7) 


If  this  differential  oquation  had  zero  on  the  right-hand  side,  it  would, 
of  course,  be  of  the  standard  Bessel's  equation  form,  of  order  iv.  As  it  is, 
it  can  be  identified  as  one  of  a  number  which  have  a  close  affinity  with 
Bessel's  form  and  which,  at  some  time,  have  received  attention  by 
mathematicians.  Equation  (7)  is  of  the  form  involving  Lommel's  functions, 
(page  40,  Vol.2  Ref. 3,  or  page  345  j  Ref.4«) 

The  complete  solution  of  oquation  (7)  will  inevitably  involve  arbitrary 
amounts  of  the  functions  J^(z),  J  ^y(z).  However,  let  us  look  at  a  particula 

solution,  (equation  (69),  page  40,  Vol.2,  Ref. 3) 


-2i\+iv 


2(i)  • 


(8) 


S.  «...  .  .  is  a  Lommel  function,  z  being  the  current  variable.  The  other 

1+2iX-iv,  iv  '•z'  *  0 

quantities,  which  are  constants,  play  a  roll  akin  to  the  order  in  normal 

Bessel's  functions. 

Lommel's  functions  (of  this  kind)  are  not  tabulated  arywhere,  as  far  as 
the  present  writer  can  discover,  so  that  one  cannot  proceed  by  consulting 
tables.  However,  some  of  their  mathematical  properties  have  been  explored, 
including  how  they  behave  for  either  small  or  large  values  of  the  argument,  z, 
which  properties  are  particularly  useful  in  the  present  context. 

By  taking  the  same  functional  form  of  the  solution  (8)  and  considering 
in  turn  its  behaviour  for  small  (z),  large  altitude  and  large  (z),  small  altitude, 
we  can  link  the  initial  re-entry  motion  at  very  high  altitude  including  a 
determination  of  the  arbitrary  pre-entry  motion,  with  the  final  stages  of  re¬ 
entry  motion  at  low  altitude  though  being  unable  to  analyse  mathematically  the 
motion  in  between.  We  will  therefore  evaluate  the  two  cases  of  low  and  high 
z  in  turn. 

3.1.1  Small  values  of  z  (large  altitude) 

Let  us  look  first  at  small  values  of  z,  corresponding  to  the  initial 
stages  of  the  re-entry  problem  when  the  altitude  is  large. 


-  7  - 


CONFIDENTIAL 


CONFIDENTIAL 


Technical  Note  No.  WE  l+O 


The  's'  function  can  then  be  expanded  in  asoending  powers  of  z. 
Equation  (8)  then  leads  to 


D  = 


-2iX+iv 


r-  z2+2iX-iv 
_2^(l +iX-iv ) (l +iX) 


4+2iX-iv  - 

z  ^ 

24  ( 1 +iX-i v ) ( 1 +iX ) ( 2+iX-iX ) ( 2+iX ) 


If  we  now  write  co  T  =  co  ,  this  reduces  to 
o  7 


a  -  ip  =  -  eT' 


2iX 


(l+iX-iv)(l+iX)  "  ( 1 +iX-iv ) ( 1 +iX) ( 2+iX-iv ) ( 2+iX )  + 


...(9) 


2iP 


2±\  iPt 

The  term  T  =  T  =  e  represents  simply  a  rotation  at  the  spin  rate,  P. 
It  means  that  equation  (9),  with  this  tern  removed,  gives  the  motion  referred 
to  missile-fixed  axes. 

Employing  a  suffix,  m,  to  indicate  reference  to  missile  axes,  we  have 


a  -  ip  = 
m  rm 


-  e 


(l+iX-iv)(l+iX)  ( 1 +iX-iv ) ( 1 +iX) ( 2+iX-iv ) ( 2+iXj  + 


(10) 


co,  which  is  the  natural  frequency  of  the  missile  in  pitch  under  aerodynamic 
forces,  gradually  increases  from  zero  during  re-entry.  Hence,  the  solution 
we  have  arbitrarily  chosen  represents  the  case  in  which  the  missile's  inertial 
axis  is  lined  up,  initially,  with  the  flight  direction,  and  commences  its 
movement  away  under  the  effect  of  the  growing,  asymmetric,  aerodynamics. 


Equation  (lO)  shows  this  in  a  typical  case,  to  be  of  the  following  form:- 


\  '  1 

\ 


-P 

-  8  - 


a 


m 

.a 

e 

i 
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The  first  terra  in  the  bracket  dominates  initially.  In  the  most  common 
case,  X  is  much  larger  than  either  1  or  v. 

In  these  circumstances, 


a  -  i6  ~  -  e 
m  rm 


1 00  \ 

V*) 


-X  +  2iX 


~  S.fS^  A  +2i\  . 

-  2{KJ  +  XJ  - 


m2  j  2i 


showing  that  the  response,  on  missile  axes,  starts  off  (for  large  X)  at  an 
angle  2/X  to  the  vector  defining  the  asymmetry. 

3.1.2  Large  values  of  z  (low  altitude) 

Now  we  turn  to  the  situation  when  z  is  large,  corresponding  to  the 
later  stages  of  re-entry,  when  the  altitude  is  much  reduced. 

We  are  dealing  of  course  with  the  same  solution,  as  defined  by 
equation  (8).  Nov/' we  wish  to  expand  in  descending  powers  of  z. 

To  do  this,  we  note  that  (equation  7*  page  4-0,  lief. 3) 


s1+2iX-iv, 


°1+2iX-iv,iv*" 

_22i\-iv  r(l+i\)r(l+i\-iv) 
sin  i-xv 


|  (l+2i\-2iv)J_iv(z) 
cos  |  (l+2i\)Jiv(z) 


(11) 


'S'  is  an  alternative  form  of  Lommel  function,  readily  expandable  in  descending 
powers.  T  is  the  gamma  function,  J  Bessel’s  function. 

For  large  z. 


^1  +2iX-iv,  iv^^ 


2iX-iv 


2iX  2i(X-v) 

2  + 
z 


±%  %v  .  ix  nv 


i"rt  %v  .  i-rc 
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When  these  expressions  are  substituted  into  (ll),  we  obtain,  after  a 
little  algebraic  simplification, 


s1+2iX-iv,iv^  ~  * 


2i\-iv 


‘  4X(X-v) 

1  +  2 

z  _J 


i2 


2iX-iv 


(2k.)1 


T"  r(l+i\)r(l+i\-iv) 


r—  iz  +  ™  k 
e  4  e 


2/  -1*  ♦  T  -  ! 

+  e  e 


•  •  •  (l2) 


Thus  from  equation  (8),  and  remembering  that 


a-  ip  =  T^Vr)  ,  z  =  “  etc.,  we  obtain 


a-  ip  ~  -  eT 


2i\ 


if  fKV/K 
2vtc  loo/  Vw 


2  2iX-iv  r  M  X 


r(l+i\)r(l+i\-iv) 


v  \  .  (  APt  2co  2 

“  V  2B  +  K  +  4, 


...  (13) 


Equation  (13)  represents  the  situation  which  ensues  at  low  altitudes. 

The  first  term  represents  a  vector  which  is  locked  to  missile  axes,  and 
rapidly  shrinking  onto  -e.  The  remainder  represents  two  vectors,  rotating 
AP 

at  rates  —  ±  00  in  relation  to  space  axes,  and  shrinking  in  inverse  proportion 
to  Vco  * 


In  general  this  is  as  far  as  the  mathematical  development  can  be  taken, 
and  tables  of  gamma  functions  can  be  consulted  for  numerical  values.  However, 
it  frequently  happens  that  X  is  much  larger  than  1  (or  v)  and  in  that  case 
approximations  to  the  gamma  functions  can  be  made. 
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The  standard  approximation  is:- 


for  large  X,  log  T(X)  ~  (X  -  -§-)log  X  -  X  +  log  2x  . 

0  6  0 


%K 


Thus, 


r(l  +  iX)  ~  .  2xX  e 


+  l 


X(log  +  | j 


I - -  |(x-v)  +  i  (x-v)  [log(x-v)  -  1}  -  p7T— r  +  r 

r(l+iX-iv)  -  \/2tc(X-v)  e  2  L  2^X_V'  4-J 


On  substituting  these  expressions  into  equation  (l3)  we  obtain,  after  a  little 
simplification, 


a  -  ip - - 


1  +  ~  (  1 
l_  a) 


ipt  1 7daTx-v ) 

+  eJ  Oi 


rHf  +  * 


/ 

-*(»-»)  +  i  (-  f 


where 

«*>  =  X  log  2X+  (X-v)  log  2(X-v) 


2X  +  V“  2X 


1  1  T. 


2 (X-v)  4 


+  7  -  (2X-v)  log 


2u> 

K 


...  (l4) 


The  equations  (14)  define  the  orientation  of  the  inertial  axis  at  low 
altitudes  for  large  X  (say  >10),  for  the  case  in  which  this  axis  had  zero 
displacement  from  the  flight  direction  initially. 

iPt 

These  expressions  refer  to  space  axes.  Deletion  of  the  terms  e  gives 

a  -  iP  ,  the  incidence  in  relation  to  missile  axes, 
m  ^m7 

3.2  General  case  including  damping  and  lift 

We  return  now  to  the  more  general  case  in  which  damping  and  lift  terms 
are  retained,  the  appropriate  differential  equation  being  equation  (5)»  ancL 
seek  a  solution  along  similar  lines  to  the  classical  mathematics  of  the  above 
special  case. 
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2a  T 

In  the  T|  term  in  that  equation,  — ~ —  will  always,  in  practice,  be  very 
4*>2T2  K 

small  in  relation  to  —  -r  —  ,  and  will  therefore  be  neglected. 

K 

We  proceed  with 


T2  *-a  +  Tin  „  f/ifa 
ai2  dT  LV 


a2T^ 

o 

K2 


-  (iv)‘ 


2  aoT 

_  e  ^  T2+i(2\-v)  JW 

K2 


(15) 


and  make  the  substitution, 


i!a  g\„ 2 

K2  '  K2  B 


a  Z 


We  then  obtain 


2  d  n  dri 
z  — i  +  z  -?*  + 
,2  dz 
dz 


z  - 


2V2  4 
a  K  z 
o 


4o2  -  2ib  ~ 
o  o  B 


-  (iv)‘ 


a0K  2 
-  z 


=  -  e 


2u)  ^  -i(2X-u)  ^  ..2  ...  AP 


2+i(2X-v)  _V“o-2ibo  By  (l6) 


to  a  sufficient  accuracy. 

If  we  let  c  =  * 


a  K 
o 


0  — Tr?r  this  equation  becomes, 

2 1  *>:-  2ibo  t 
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2  dj 

5  7*2 

dz 


+  z 


drj 

dz 


2,24  /.  n2 

z  -  4c  z  -  (iv) 


T]  =  -  e 


2b) 

< 

K 


-i(2X-v) 


2+i(2X-v)  cz' 
z  '  '  e 


...  (17) 


This  is  the  counterpart  to  equation  (7)  in  Section  3*1 «  The  term 

2  4 

4c  z  in  the  coefficient  of  r]  is  negligible,  in  a  typical  re-entry  case,  in 

comparison  with  z  .  It  is  best  retained,  however,  since  it  assists  the 
mathematical  development. 

3.2.1  Ascending  powers  of  z 

Following  the  lines  of  Section  3*1 >  we  seek  first  a  solution  in 
ascending  powers  of  z  corresponding  to  Lommel's  function,  s.  Denoting  the 
corresponding  function  by  "s,  and  writing  |i  =  1+i(2\-v),  we  can  devise 


s  =  e 


cz 


u+1  u+3  u+5 

£t .  Z  4*  Q,  7  Z  +Q.J-Z  •  •  •  • 

1  3  5 


in  which 


(ji+l)“-  (iv)2 


1  +  4c(u+2) 


[(|i+l)2-  (iv)2][fi+3)2-  (iv)2] 


[1  +  4c(u+2)][l  +  4c(u+4)] 


5  [(t.+l)2-(iv)2][(M+3)2-(iv)Z][(|Jt5)2-(iv)Z] 


;  etc. 


,2u 

1  ?1+2i%-iv,iv<l)  (l8> 


corresponding  to  the  case  in  which  the  inertial  axis  is  lined  up  initially 
with  the  flight  direction. 

3. 2. 2  Descending  powers  of  z 

Here  we  seek  a  form  of  solution  to  equation  (l7)  corresponding  to 
Lommel's  function  S.  Denoting  this  by  S,  we  can  obtain 
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in  which 


ao  “  1+4cji  '  a2 


[(u-l)2-  (iy)2  3 

[l+4cnj[l+4c(|i-2)  j 


Thus, 


-  e 


/2w 


-i(2\-v) 


1+2iX-iv,iv 


(z) 


providing  a  convenient  form  of  solution  for  large  z. 

3.2.3  Complementary  function 

We  require  a  suitable  form  for  the  complementary  function  of  equation  (l7)» 
i.e.  the  solution  when  the  right-hand  side  is  set  to  zero,  and  explore  the 
form 


f] 


2 

oz  r  r  r+2 

e  |^aoz  +a2z 


+  a.  z 
4 


r+4 


+ 


By  substitution  in  (17),  it  is  found  that 

or  r  =  -iv 

arbitrary 

a  [l  +  4c(-iv+1 )] 

2  (-iv+2)2-  (-iv)2 

aQ[l  +  4c(iv+1 )][l  +  4c(iv+3)]  aQ[l  +  4c(-iv+1 )][l  +  4c(-iv+3)] 

^  [(iv+2)2-  (iv)2][(iv+4)2-  (iv)2]  ^  [(-iv+2)2-(-iv)2][(-iv+4)2-(-iv)2] 

etc.  etc. 


r  =  +iv 
aQ  arbitrary 

aQ[l  +  4c(iv+  l)] 
2  (iv+2)2-  (iv)2 
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Thus,  the  two  solutions  of  this  form  can  be  written 


A  eCz2  -  t1  +  4c(iv+1  )]z2  _ 

'  22(iv+1 )  jl  24(iv+l)(iv+2)|2 


(19) 


A,  eCZ^  z"iv f  1  -  L1  .tAc(-ivtl),]^  +  [1  +  4c(-iv+l)][l  +  4c(-iv+3)jz4 


2(-iv+1 )  Li.  24(-iv+1 )(-iv+2)  |2 


J 


(20) 


where  A^  and  A^  are  arbitrary  constants. 


As  in  the  standard  Bessel  case  (corresponding  to  c  =  0),  we  take 


A 


1  „iv 


Ar 


2  r(iv+l) 


2_lv  r(-iv+l) 


.  ( cf .  kef. 4, p» 40) 


Q  I 

Now,  in  re-entry  applications  the  quantity  c  is  quite  small,  and  also  or z  is 

2 

always  quite  small  in  relation  to  z  .  Thus  for  the  range  of  z  in  which  we  are 
interested,  equation  (17)  with  the  right-hand  side  set  to  zero  is  sensibly 
equivalent  to  the  standard  Bessel  form  of  order  iv,  and  the  complementary 
function  solutions  found  can  be  identified  for  our  purpose  as  J^(z)  and 

J  ^  (z)  respectively. 


3.2.4  Complete  solution 

The  solution  in  which  we  are  interested  is  as  given  by  equation  (l8) 
since  this  corresponds  to  zero  initial  displacement.  We  require  now  to 
express  s  in  terms  of  S  and  the  complementary  function  <L^  and  J  .  .  This 

will  be  the  counterpart  in  the  more  general  case  to  equation  (ll)  in  the 
special  case,  and  will  enable  an  asymptotic  value  to  the  solution  to  be  found. 
This  in  turn  will  be  used  to  explore  what  happens  to  the  motion  in  the  later 
stages  of  re-entry. 

It  is  a  matter  of  algebra  to  show  that  there  is  an  analogous  relation 
to  equation  (ll),  viz:- 
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s1+2iX-iv,iv^  =  ^1+2iX-iv,iv^ 


22i\-iv  r(uiX)r(l+iX-iv) 


sin  ixv 


* 


J_  lx  ir 

^8c+  2  iX"  2  j 


'^UiX27  c°.i(^21X-2i,)  J.iv(2) 


r  i-L+i+iy 

\8c  2  2 


(8c) 


1+iX-iv  008  2  (U21X)  Jiv(z) 


(21) 


where  expression  (l9)  and  (20)  are  used  for  J^y(z),  J_iy(z)  in  this  case 


-iv' 


As  in  Section  3*1*2,  in  order  to  proceed  further  with  the  mathematical 
development  it  is  necessary  to  approximate  to  the  gamma  functions.  This  is 
quite  feasible,  since  in  practice,  ’o’  is  quite  small  and  l/o  quite  large. 

Using  the  standard  approximation  for  the  T-function  of  a  large  quantity, 
it  is  straightforward  to  show  that  for  Y  large,  and  much  larger  than  z 


By  making  use  of  this  formula,  it  can  be  demonstrated  that 


r  ( J 

\8 


8c 


r 


L  +  2 

•8o  2 


1 

(8c)1+iX 


j  1  1  iv ) 

1  V>8o  *  2  *  2  J  1 

„  / 1  3  iv\  /a  \  1+iX-iv 

r  (to  *  2  +  iX  "  —)  (8o) 

a  K 

[X(X-v)  +  i(-2X+v)] 


Also  from  Section  3*2.2,  we  have,  for  large  z, 


1 +2iX-iv,iv 


(.)  e0"2  z^-1  1  -  TIie=ljfLJi>dfL  _L  +  o  (X\ 

'  ’  ~  |_1+4c|i  [  1  +4c|a J  [  1  +4c ( (a- 2 )  J  2  \2V_ 


2 

cz 
e  z 


i(2X-v)  [~  _  .  tli.  ( X  -  ^  +  &i±=vl 

L1  a.2  2u)2  v  y  z2  . 
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Hence,  [c.f.  equation  (l2).  Section  3*  1«2] 


-  /  \  cz  i(2\-v) 

s.  .  .  U)  ~  e  z  ' 

1+2i\-iv,ivN 


i —  a  K  ia  K  /  \  . « /.  \ 

A  _  -2-  _  — £L_  A  _  v\ 

1  2  ,2  r  2)  2 

L  2co  2co  \  /  z  _l 

cj  o 


,,l(2X-v)  tX(X-v)+i(-2Uy)] 

+  — - r—  r(l+i\)r(l+i\-iv)  e  0 


(2n:z): 


*(h)  -*(H) 

6  ex  '  4.  0  ^0  N  ' 


]• 
...  (22) 


Recalling  that 


T|  = 


e 


(2X-v) 


s 


> 


a  -  ip 


H  > 


v2  2 

K  z 


T 


2 


and  using  approximations  to  the  functions  F(l+i\),  r(l+i\-iv)  as  in  Section 
3.1.2  we  can  find  the  asympotic  value  of  a-ip  corresponding  to  late  stages  of 
re-entry. 

This  comes  out  to  be 


a-ip - - 


K 


a 


2  2 

0)  CO 

o  o 


1  +  2  --3+-75V  X-S  -  i  X  -  £ 


(0 


p 


CO 


1  -  ~ 


iPt 


a  KX(X-v) 


+  ,  2ja(X-v), 

J  CO 


2co 


o  2K 

e 


a  bv  .  /  2co 

r%  +  1  It  + 


K 


+  e 
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where 


<t>  =  X  log  2X  +  (X-v)  log  2(X-v)  -  2X  +  v  -  ~ 


-  (2X-v)  log 


a  K 

“*5  (~2X+v) 


b  Kv 

+  "V 

8(0 

0 


(24) 


This  is  the  mathematical  solution  of  the  problem,  subject  to  certain 
restrictions  on  the  size  of  some  of  the  quantities  involved  (see  Section  4)« 
It  is  the  more  general  form  of  equation  (l4) >  the  effect  of  damping  and  lift 
introducing  the  extra  terms  in  a,  b,  aQ,  bQ. 

4  PRACTICAL  IMPLICATIONS 


We  have  found  [equations  (23),  (24)3  the  angular  motion  of  an  asymmetric 
spinning  body,  occurring  at  low  altitudes,  as  a  consequence  of  its  inherent 
asymmetry.  As  appears  later  in  this  paragraph,  this  is  the  most  suitable  form 
for  an  appraisal  of  the  practical  implications  of  asymmetry.  For  equations 
(23)  and  (24)  to  be  valid,  two  conditions  must  be  satisfied: - 


a  b 

(a)  damping  and  lift  must  be  relatively  small,  i.e,  —z  ,  must  be 

co 


small,  in  comparison  with  unity, 


00 


must  be  large  in  comparison  with  unity. 


For  all  typical  re-entry  cases,  the  former  condition,  (a),  is  well  satisfied. 
For  condition  (b),  it  is  sufficient  if  X  is  greater  than  about  10,  and  this 
is  frequently  the  case  in  re-entry  applications.  If  it  should  happen  that  X 
is  much  less  than  this,  it  would  be  necessary  to  go  back  to  equation  (22)  and 
evaluate  the  T-functions  from  suitable  tables. 

The  effect  of  asymmetry,  as  per  equation  (23),  consists,  in  the  later 
stages  of  re-entry,  of  two  parts 

(l)  a  steady  displacement  of  the  inertial  axis  of  amount 


00  0 


in  relation  to  missile  axes. 

This  is  approximately  equal  and  opposite  to  the  aerodynamic  asymmetry, 

(2)  a  motion  consisting  of  two  rotating  vectors  of  exactly  the  same 
type  as  the  normal  symmetrical  missile  undergoes. 
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If  e  is  small,  as  is  supposed,  then  constituent  (l)  is  also  small.  Interest, 
in  practice,  centres  on  constituent  (2).  Of  this  the  component  vector  which 
rotates  in  the  same  sense  as  the  missile  spins  will  be  much  larger  than  the 

counter  rotating  vector,  since  the  latter  involves  the  factor  e  7':(2X“V)# 

The  feature  of  principal  interest,  therefore,  is  the  vector 


aKXiVvie 

0) 


aQKX(X-v) 

“£1?— 


a  bv  .  ( 2o> 

+  —  +  i  [  y  +  <P 


2K  2u) 


Now,  it  is  shown  in  Ref.1,  that  if  a  symmetrical  missile  has  a  simple 
precessional  motion  of  magnitude  S^,  prior  to  re-entry,  then  its  motion  during 

re-entry  will  be  given  by 


a  bv 

q  j  AP  2K  +  2u)  +  1 
E  J  2Bw  6 


2o) 

K 


+  $ 


These  two  expressions  are  of  the  same  foim.  Hence,  by  equating  'them, 
we  oan  express  the  effeot  of  asymmetry  in  terms  of  an  equivalent  motion,  prior 
to  re-entry  for  a  symmetrical  missile 

a  KX(X-v) 


E 


2oo 


=  e  e 


>1 


xia(x-v) 


Therefore 


SE 


a  KX(X-v) 

Q 

2o)2 


=  e 


27i;X(x-v  \ 

V 


(25) 


For  this  relationship  to  be  valid,  the  further  condition  that  v  must  not 
be  less  than  about  unity  must  be  satisfied.  This  is  usually  so  in  practice. 

Rewritten  in  terms  of  basic  design  parameters  equation  (25)  becomes 


=  e 


#>2(cr2C  -  C  )  (1  -  A/B) 

a  q  _ 

2KV„  C 
E  m 
a 


2x 


P  /B 
K  \  A  “ 


(26) 
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where 


O' 


P 

A,B 

VE 

K 


aerodynamic  derivatives  of  lift,  pitch  damping,  and  static 
pitching  moment 

reference  length 

non-dimensional  radius  of  gyration  in  pitch,  i.e.  ratio  of 
pitch  radius  of  gyration  to  reference  length 

spin  rate 

moments  of  inertia  of  re-entry  missile  in  roll,  pitch 
re-entry  velocity 

Kt 

given  by  air  density  variation,  p  =  p  e 
asymmetry  in  missile 

equivalent  re-entry  perturbation  for  symmetrical  missile 


Equation  (26)  is  the  principal  practical  result  of  the  present  enquiry. 

It  defines  the  situation  in  which  the  missile  finds  itself  after  passing 
through  resonance,  and  expresses  it  in  terms  of  the  equivalent  re-entry  dis¬ 
turbance  for  a  symmetrical  body. 

It  is  interesting  to  note  that  the  presence  of  damping  (and  lift)  in 
the  re-entry  missile  increases  the  magnification  factor  (as  defined).  For 
positive  lift,  pitch  damping,  and  static  pitching  moment,  is  positive, 

a 

and  Cm  negative,  so  that  the  exponential  term  in  (26)  consists  of  'e'  to 
q  a 

a  positive  power,  and  this  is,  of  course,  greater  than  1. 

It  will  usually  be  the  case,  however,  that  the  exponent  of  'e'  will  be 
fairly  small  -  a  value  of  0*1  to  0*2  would  be  typical  of  ourrent  design 
practice  -  so  that  the  contribution  of  the  exponential  term  in  (26)  will 
frequently  not  be  very  important. 

The  main  effect  will  generally  come  from  the  remaining  term  in  equation 

-  l\  .  In  current  designs  this  tends  to  be  of  the  order 

of  30  to  40.  This  means,  therefore,  that  if  such  a  re-entry  missile  is  asym- 
metrio  to  the  extent  of  1°  angular  difference  between  the  aerodynamio  and 
inertial  axes  and  it  has  no  precession,  the  effeot  on  the  incidence  in  the 
later  re-entry  stages  is  the  same  as  if  the  missile  were  symmetric  but  had  a 
re-entry  incidence  in  the  form  of  a  precession  of  30°  to  40°  semi-angle 
round  the  flight  path. 


(26),  namely 


*1(1 


-  20  - 


CONFIDENTIAL 


CONFIDENTIAL 


Technical  Note  No.  WE  40 


5  NUMERICAL  CHECK  ON  THEORETICAL  RESULTS 


It  was  considered  desirable  to  check  the  validity  of  the  foregoing 
theoretical  analysis  for  one  or  two  particular  cases.  For  this  purpose,  a 
digital  calculating  machine  was  programmed  to  solv$  the  basic  differential 
equation,  (4).  There  is  one  minor  difficulty  which  had  to  be  overcome  -  not 
in  "the  programming  as  such,  but  in  feeding  in  appropriate  initial  conditions* 
The  mathematical  solution  corresponds  to  initial  conditions 


a  -  ip 


a  -  ip 


0  at  t  =  -00  , 


To  overcome  this,  values  for  a-  ip,  a  -  ip  were  calculated  for  a  later  time, 
appropriate  to  the  very  early  stages  of  re-entry,  from  the  expression  for 
s  (or  s),  and  the  digital  calculation  proceeded  from  there. 


Fig.1  shows  the  digital  machine  output  for  three  particular  cases,  with 
progressively  increasing  damping.  For  each  of  them 


P  =  20,  K  =  1,  A/B  =  0  . 


The  damping  corresponds  to 


a 

b 

a 

b 

a 

0 

_  _o 

s  0  .  -2  s 

_ 0  _  i_ 

0 

2 

2 

’  2 

n  ~  A  9 

2 

0) 

0) 

0) 

0) 

0) 

0 

0 

0 

0 

0 

1 

*  ?  > 


the  lift  coefficient  being  zero  in  each  case. 


The  first  one  is  the  case  of  zero  damping.  The  other  two  are  equivalent 
to  rather  more  damping  than  current  re-entry  designs  normally  possess,  but  it 
was  desired  to  accentuate  the  effect. 


The  mathematical  solution,  as  per  equations  (23),  (24),  was  checked 
against  the  digital  machine  results.  The  agreement  is  virtually  perfect  for 
w/P  greater  than  1*5,  and  is  still  very  good  for  co/P  as  low  as  1  *3* 

Other  cases  were  subjected  to  the  same  check,  with  different  values  of 
A/B  and  damping  and  lift.  They  present  the  same  general  appearance,  and 
confirmation  of  the  mathematical  theory. 

6  CONCLUSION 


A  mathematical  solution,  in  convenient  form,  has  been  derived  for  the 
effect  on  the  angular  disturbance  of  a  spinning  re-entry  missile,  which  will 
result  from  its  having  a  small  amount  of  asymmetry. 
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The  solution  is  readily  applicable  in  practice,  since  in  the  post¬ 
resonance  region  on  which  interest  usually  centres,  the  effect  is  simply 
related  to  the  case  of  the  corresponding  symmetrical  missile. 


A,B 


a,b 

c 

K 

l 

m 

P 


q,r 

S 


NOTATION 

moments  of  inertia  in  roll  and  pitch  (or  yaw) 

parameters  defining  damping  and  lift,  Section  2 

parameter  relating,  inter  alia,  damping,  lift  and  frequency 

Kt 

parameter  defining  variation  of  air  density  with  time,  p  =p  e 

o 

aerodynamic  reference  length 
mass 

spin  rate 

pitch,  yet w  angular  rates 
aerodynamic  reference  area 


a  a  q 


a,p 

Y>6 

e 

X 


Lommel  functions  (usually  with  two  suffixes) 

Kt 

2 

e 

time 

re-entry  velocity 
aerodynamic  derivatives 

missile  incidence  components 
trajectory  deflection  angular  components 

missile  asymmetry 

P 

K 

1  +  i(2X-v) 


v 


AP 

BK 
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NOTATION  (CONTP) 


p  air  density 

co  parameter  defining  frequency,  Section  2 

Suffix 

o  condition  at  t  =  0 
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